clear;
load('datamarketG2')


dt=0.01;
T1=30;
t1=0:dt:T1;
nt1=length(t1);



gnmk2=phiXmk2(1:flagstarmk2)-phiImk2(1:flagstarmk2);





nmk2(1)=ss(1);
for i=2:nt1
nmk2(i)=nmk2(i-1)+gnmk2(floor((nmk2(i-1)-ss(1))/ds)+1)*dt;
end






 
 plot(t1, nmk2,'r--','LineWidth',1.5)
 hold on;

 
 axis('square')
 
  xlabel('$t$','interpreter','latex','FontSize',12)

  
  
  

title('$n_t$','interpreter','latex','FontSize',12)
 
 axis('square')

 reportnew2;
 
 print -depsc figure_base_transition_new13_power_zeta02_new_stateG55